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Abstract. We present accretion rates obtained from three-dimensional self-gravitating radiation 
hydrodynamical models of giant planet growth. We investigate the dependence of accretion rates 
upon grain opacity and core/protoplanet mass. The accretion rates found for low mass cores are 
inline with the results of previous one-dimensional models that include radiative transfer. 
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INTRODUCTION 

The core accretion model of planet growth has been devised to explain the formation of 
giant planets such as Jupiter and Saturn [1, 2]. Through the coagulation of solid grains, 
planetesimals are built which then, through a so far uncertain process, bind to form large 
km scale bodies. When such bodies collide they shatter, but more often than not are held 
together by their self-gravity [3], allowing them to grow through collisions to reach many 
Earth masses. These large solid masses exert a considerable gravitational tug on the solar 
nebula in which they are immersed, causing them to slowly gather a gaseous atmosphere. 
If the solid cores are sufficiently massive (> 10 MEarth) they may trigger runaway gas 
accretion. This occurs when the accreted gaseous atmosphere obtains a mass similar to 
the solid core, causing it to collapse from its hydrostatic state. The collapse leads to a 
rapid in fall of material from the surrounding nebula to replace the condensing gas. It is 
believed that such a process allows the accretion of hundreds of Earth masses of gas. 

Three-dimensional hydrodynamical calculations examining the interaction of a pro- 
toplanet with a disc have usually assumed a locally isothermal equation of state, and/or 
do not resolve the gas flow far inside the planet's Hill radius [4-13]. Recent hydrody- 
namical calculations have also considered different equations of state and/or radiation 
transport [10, 12, 14-16]. Fouchet & Mayer [17] include both radiative transfer and self- 
gravity, but once again the resolution of gas flow near to the planetary core is limited. 
In this work we present results from three-dimensional self-gravitating radiative transfer 
hydrodynamical calculations that are able to resolve gas flow down to realistic core radii. 
We vary the core/protoplanet masses and the grain opacity of the nebula to examine the 
effects that these properties have upon the gas accretion rates of the growing planets. 

IMODEL SETUP 

We model 10 - 333 Mearth solid cores/pro toplanets embedded in a circumstellar disc 
with a surface density of 75 g cm~^ orbiting a IM0 star at 5.2 AU. The initial surface 



density has a E r^^l^ profile, whilst the initial temperature profile of the disc is 
r oc r 1, yielding a constant disc scaleheight with radius oi H /r — 0.05. We vary the 
grain opacity of the disc from the interstellar grain opacity (IGO) down to 0.1% of this 
value to examine the effects this has upon the disc cooling. 

To achieve high resolution near a protoplanet at reasonable computational cost we 
model only a small section of the protoplanetary disc centred on the planetary core. Our 
section size is 5.2± 0.78 AU, and = ±0.15 radians; for details see [18]. 

RESULTS & DISCUSSION 

The accretion rates for the different mass cores/protoplanets that we modelled under 
differing opacities are shown in the left panel of figure 1. The trend in accretion rate 
against core/protoplanet mass is as would be expected, increasing as the central bodies 
mass increases. There is a reduction in accretion rates for the highest mass protoplanets 
caused by the evacuation of a large gap in the encompassing circumstellar disc by 
these protoplanets. Such a gap reduces the material available in the feeding zone of 
the protoplanet. The accretion rates also show an increasing trend with reduced grain 
opacities for the low mass cores. Reduced grain opacities allow the infalling material 
to cool more readily via radiation, and so the thermal pressure supporting it against 
capture is reduced. The only exception to this is at the lowest grain opacity of 0.1% 
IGO for the lowest mass core where the radiation and gas fields interact so little that the 
two fields decouple, preventing effective gas cooling via radiation. For the high mass 
protoplanets the grain opacity has very little impact on the accretion rates as the process 
is gravitationally dominated. 

The right hand panel of figure 1 compares the accretion on to a lOMEarth solid core in 
our models with the accretion onto 5 and 15 MEanh cores in the one-dimensional models 
of Papaloizou & Nelson [19]. It can be seen that the results show good agreement, sug- 
gesting that for at least the phase of accretion preceding runaway growth, the transition 
from one-dimensional models to three-dimensional models does not have a significant 
impact on the accretion rates. 

Figure 2 illustrates the differences between radiation hydrodynamical and isothermal 
calculations. Features in the vicinity of the planet that are well established under isother- 
mal conditions, such as spiral shocks and circumplanetary discs, are weakened in the 
radiative transfer calculations. The thermal pressure acts to smear out these features. It 
can also be seen in figure 2 that reducing the grain opacity leads to a model that is more 
similar to the isothermal case. The increased cooling seen at lower opacities reduces the 
thermal pressure, thus mitigating its impact on the shocks and circumplanetary discs. 
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FIGURE 1. a. Protoplanet accretion rates. The asterisks mark the accretion rates for our SPH calcula- 
tions using a locally-isothermal equation of state, essentially providing upper limits to the accretion rates 
obtainable with various protoplanet masses. The diamonds mark the results of [9], and the plus signs 
those of [6], connected with solid lines. These calculations were also locally-isothermal but were global 
disc simulations performed using the ZEUS code. The SPH and ZEUS accretion rates are in reasonable 
agreement. The accretion rates from our self-gravitating radiation hydrodynamical SPH calculations using 
core radii of 1% of the protoplanet' s Hill radius are given using line types that denote the different grain 
opacities. Results are shown using standard interstellar grain opacities (IGO) (dotted), 10% IGO (short- 
dashed), 1 % IGO (dot-dash), and 0.1% IGO (long dashed). The inclusion of radiative transfer substantially 
lowers the accretion rates of low-mass protoplanets, but Jupiter-mass protoplanets have similar accretion 
rates to the locally-isothermal result, regardless of the grain opacity. The analytic approximation of [20] 
is shown by the solid curved Une. b. Accretion rates versus the gas mass accreted by planetary cores for 
1% (thin lines) and standard (thick lines) grain opacities. We plot the results of [19] for 5 ME^nh (dashed 
Unes), and 15 M^Eanh (dot-dashed lines) cores. The soUd Unes show the extrapolations of our 10 M-Earth 
core results. Our 10 Mgarth core accretion rates lie between the accretion rates obtained by Papaloizou 
and Nelson for their 5 MEanh and 15 ^Earth cores. The effect of the reduced opacities is also consistent 
with Papaloizou and Nelson's results. 
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FIGURE 2. Surface density plots for locally-isothermal calculations (left), and self-gravitating radiation 
hydrodynamical calculations using 1% IGO (centre), and standard IGO (right) calculations with our 
standard protoplanetary disc surface density. From top to bottom the protoplanet masses are 22, 33, 100, 
166 and 333 MEaith respectively. The radiation hydrodynamical calculations use protoplanet radii of 1% 
of the Hill radii, while the locally-isothermal calculations use 5% of the Hill radii. Note that with radiative 
transfer and standard opacities, the spiral shocks in the protoplanetary disc are much weaker than using 
a locally-isothermal equation of state, while using radiative transfer with reduced grain opacities results 
in intermediate solutions because the discs are less optically thick and are able to radiate more effectively 
than with standard opacities. 



